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Resumen 


Se presenta un modelo numérico para riego por melgas basado en el acoplamiento 
interno de las ecuaciones de Saint-Venant y Richards, el cual permite describir las 
fases de avance, almacenamiento, consumo y recesión del riego por melgas. La 
solución de las ecuaciones de Saint-Venant se aproxima utilizando un esquema 
Lagrangiano en diferencias finitas, mientras que para la ecuación de Richards se 
utilizan elementos finitos para la integración en el espacio y diferencias finitas 
implícitas para la integración en el tiempo. Los resultados obtenidos permiten 
concluir que el modelo puede ser utilizado con fines de diseño de riego por 


melgas. 


Palabras clave: fase de avance, fase de almacenamiento, fase de consumo, fase de 


recesión, acoplamiento numérico interno. 


Introducción 


La descripción del flujo del agua en el riego por 
melgas requiere tener en cuenta las ecuaciones 
de Saint-Venant para describir el flujo del agua 
sobre la superficie del suelo y la ecuación de 
Richards para modelar el flujo del agua en el 
suelo, de tal manera que su acoplamiento per- 
mita conocer las formas de los perfiles de flujo 
superficial y la distribución de los potenciales 
de presión en el suelo según avanza el riego. 
El tirante proporcionado por las ecuaciones de 
Saint-Venant se utiliza para definir una condición 
de frontera tipo Dirichlet para la ecuación de 
Richards, cuya solución numérica permite el 
cálculo de la lámina infiltrada, necesaria para 
la solución de las ecuaciones de Saint-Venant, 
teniéndose, por lo tanto, un procedimiento 
iterativo. 

Los principales trabajos de modelación 
del riego por gravedad pueden diferenciarse 


según el tipo de ecuaciones que emplean 
para representar el fenómeno. Respecto a la 
modelación del flujo de agua a superficie libre 
sobre el suelo, se emplean básicamente cuatro 
aproximaciones: modelo hidrológico o de balan- 
ce de volumen, onda cinemática, onda difusiva 
y ecuaciones de Saint-Venant completas. Para la 
representación del flujo del agua en el suelo se 
utilizan principalmente tres aproximaciones: la 
ley de infiltración empírica de Kostiakov (1932), 
leyes de infiltración con base físico-matemática 
e hipótesis simplificadoras: ecuación de Green y 
Ampt (1911), ecuación de Philip (1957) y ecuación 
de Parlange et al. (1982), y ley de infiltración 
conforme a la ecuación de Richards (1931). 

En la literatura relacionada con el tema de 
este trabajo se reportan diversos modelos para 
simular el flujo del agua en el riego por melgas, 
los cuales difieren por la forma en que se 
aproxima la descripción tanto del flujo del agua 
a superficie libre como del flujo del agua en el 
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suelo, en este contexto pueden citarse trabajos 
en los que el comportamiento del medio po- 
roso se modela mediante leyes de infiltración 
con formas potenciales, como la ecuación de 
Kostiakov-Lewis (Sakkas y Strelkoff, 1974; 
Strelkoff y Katopodes, 1977; Katopodes y 
Strelkoff, 1977; Abiola y Katopodes, 1988; 
Catalán et al., 1988; Wallender y Rayej, 1990; 
Bautista y Wallender, 1992; Playán et al., 1994; 
1994; Vivekanad y Bhallamundi, 
1995), o bien utilizado la ecuación de Parlange 
(Schmitz et al., 1985; Vivekanad y Bhallamundi, 
1995) o la ecuación de Philip (Philip y Farrell, 
1964; Or y Walker, 1996). Recientemente se han 
desarrollado modelos que utilizan ecuaciones 


García, 


diferenciales bidimensionales para describir 
el flujo del agua en el suelo (Bradford y 
Katopodes, 1998; Saucedo et al., 2005; Banti 
et al., 2007; Bautista et al., 2010). En relación 
con los dos últimos trabajos citados, los más 
recientes, Banti et al. (2007), proponen un 
modelo basado en el uso de las ecuaciones 
de Saint-Venant y Richards, pero difieren 
significativamente al desarrollado en este 
trabajo en la forma de resolver numéricamente 
las ecuaciones de Saint-Venant, ya que en el 
artículo citado se emplea un esquema numérico 
tipo MacCormack, cuya aplicación requiere 
del uso de la hipótesis de un tirante pequeño 
presente sobre toda la melga desde el inicio 
del riego, lo cual no guarda concordancia con 
la física del riego por melgas; en el presente 
trabajo se prescinde de dicha hipótesis, al 
emplear un método euleriano-lagrangiano 
basado en Strelkoff y Katopodes (1977), el cual 
permite seguir detalladamente la evolución 
del frente de avance. El trabajo de Bautista et 
al. (2010) presenta un análisis de la naturaleza 
del acoplamiento de las ecuaciones empleadas 
para modelar el riego por melgas, y propone el 
empleo de un acoplamiento externo, diferente 
al acoplamientointernoutilizado en este trabajo; 
según los autores citados, el acoplamiento 
externo reduce el tiempo de cómputo y mejora 
la convergencia. En comparación con el trabajo 
presentado, puede decirse que con el empleo 
de nuevas computadoras es posible reducir 


el tiempo de cómputo, y que la mejora en la 
convergencia puede alcanzarse reduciendo 
los pasos de tiempo y espacio en el modelo 
presentado en este trabajo. 

En el riego por melgas pueden distinguirse 
cuatro fases: avance, almacenamiento, consu- 
mo y recesión. La fase de avance inicia 
con la introducción del agua en la melga y 
termina cuando el frente de onda alcanza el 
extremo aguas abajo de la melga. La fase de 
almacenamiento comienza cuando el frente de 
onda llega al extremo aguas abajo de la melga 
y culmina cuando se deja de aplicar el gasto de 
aporte en la cabecera. La fase de consumo inicia 
cuando se deja de aplicar el gasto de aporte en 
la cabecera de la melga y culmina cuando se 
tiene un tirante nulo en dicha posición. La fase 
de recesión comienza a partir de que se tiene un 
tirante de agua nulo en la cabecera de la melga 
y termina cuando se alcanza un tirante de agua 
nulo en el extremo aguas abajo de la melga. 
El objetivo de este trabajo es desarrollar un 
modelo basado en el acoplamiento numérico 
interno de las ecuaciones de Saint-Venant y 
Richards que permita describir todas las fases 
del riego por melgas. 


Flujo del agua sobre la superficie del 
suelo 


El flujo del agua con una superficie libre 
es modelado con las ecuaciones de Saint- 
Venant, que resultan de la aplicación de las 
leyes de conservación de masa y cantidad 
de movimiento. En una melga, la relación 
entre su ancho y el tirante de agua permite 
considerar las ecuaciones correspondientes al 
escurrimiento sobre una superficie de ancho 
infinito (Woolhiser, 1975). La ecuación de 
continuidad se escribe como: 


A A 


= 1 
ot ox ot (1) 


La ecuación de momentum se escribe en la 
forma recomendada por Saucedo et al. (2005): 
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hot n? ox nh? lox 
(2) 
q ol 
——=10) 
+ 


donde q(x,t) = Ul(x,t)h(x,t) es el gasto por 
unidad de ancho de melga [L?T"]; x, la 
coordenada espacial en la dirección principal 
del movimiento del agua en la melga [L]; t, el 
tiempo [T]; U, la velocidad media; h, el tirante 
de agua [L]; J, la pendiente topográfica de la 
melga [LL”]; J, la pendiente de fricción [LL”*); 
V, = 01/0t, el flujo de infiltración [LT”'], es decir, 
el volumen de agua infiltrado en la unidad de 
tiempo por unidad de ancho y por unidad de 
longitud de la melga; l, la lámina infiltrada [L]; g, 
la aceleración gravitacional [LT”]; el parámetro 
adimensional f = 1 -— a, siendo a. = 1 - U, /U, 
donde UI, es la proyección en la dirección del 
movimiento de la velocidad de salida de la 
masa de agua debido a la infiltración. 

La relación entre las variables hidráulicas q 
y h con la pendiente de fricción, denominada 
ley de resistencia hidráulica, es adoptada de 
acuerdo con Fuentes et al. (2004), haciendo uso 
de una ley potencial de resistencia: 


d 
hs 
=kv 


donde v es el coeficiente de viscosidad cine- 
mática del agua [L?T"] y k es un factor 
adimensional. 


(3) 


Para establecer las condiciones inicial y de 
frontera que deben sujetar a las ecuaciones de 
Saint-Venant para modelar las cuatro fases del 
riego por melgas, se tiene en cuenta una melga 
cerrada, a fin de evitar la pérdida de masa fuera 
del dominio de riego. Para la fase de avance se 
tiene: 

H(x,0)=0 


q(x,0)=0 y (4) 


q0,t)=3,  qlept)=0 y  h(xpt)=0 (5) 


donde x(t) es la posición del frente de onda 
para el tiempo t y q, es el gasto de aporte 
en la entrada de la melga. Para la fase de 
almacenamiento: 
q(0,t) =9,, q(L,t)=0 (6) 
donde L es la longitud de la melga y q, es el 
gasto de aporte en la cabecera de la melga. Para 
la fase de consumo: 
q(0,)=0, — (L,A=0 (7) 
Para cerrar el sistema es necesario conocer 
la forma en que evoluciona en el tiempo la 
lámina infiltrada en toda posición sobre la 
melga, es decir, la ley de infiltración, la cual 
es proporcionada con base en la solución de la 
ecuación de Richards (1931). 


Flujo del agua en el suelo 


La descripción del proceso de infiltración del 
agua en el suelo puede efectuarse teniendo 
como base la ecuación de Richards (1931), 
que resulta de la combinación de la ecuación 
de continuidad con el campo de velocidades 
calculado conforme a la ley de Darcy y que en 
su forma tridimensional, sin tener en cuenta la 
extracción de agua por las plantas, se escribe 
como: 


Cy) a =V-[K(wvy]- EY) 
5 y dz 


donde y es el potencial de presión del agua 
en el suelo, expresado como la altura de una 
columna equivalente de agua [L] (positivo en la 
zona saturada y negativo en la zona no saturada 
del suelo); C(w) = d8/dy es denominada la 
capacidad específica de humedad del suelo; 
9 = 0 (y) es el volumen de agua por unidad 
de volumen de suelo o contenido volumétrico 
de agua [L*L*] y es una función de y conocida 
como curva característica de humedad o 
curva de retención del agua; K = K(y) es la 
conductividad hidráulica [LT“], que en un 
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suelo parcialmente saturado es una función del 
potencial de presión; el potencial gravitacional 
es asimilado a la coordenada espacial z 
orientada positivamente hacia abajo [Ll]; 
V = (9/0x,9/0y,0/0z) es el operador gradiente; 
x, y son las otras dos coordenadas espaciales 
[L], y t es el tiempo [T!]. 

El riego es un fenómeno que se efectúa en 
tres dimensiones, lo cual hace necesario que el 
flujo del agua en el suelo sea descrito utilizando 
la ecuación (8), sin embargo, es conveniente 
aceptar la hipótesis de que el fenómeno se 
efectúa en planos paralelos al desarrollo de 
la melga; entonces es posible utilizar la forma 
bidimensional de la ecuación de Richards: 


dy _9d oy| 0 dy 
me? ES y . 2 EN o J] 0) 


la cual debe resolverse sobre el dominio de 


solución que se muestra en la figura 1. 

Como condición inicial para la solución de 
la ecuación de Richards bidimensional se debe 
especificar la distribución de las presiones en 
el espacio: 

y =v, (%2) (10) 

Las condiciones de frontera correspondien- 
tes pueden considerarse como sigue: Ax; 


frontera tipo Dirichlet con potencial prescrito 
conforme a la solución de las ecuaciones de 
Saint-Venant, x,B, BC y DA fronteras tipo 
Neumann con flujo nulo, CD frontera bajo 
gradiente unitario: 


y=), x€ AXy, z=0, t>0 (11) 
x= o, xexB, z=0, t>0 (12) 
Z 
dy — 
Kw)> =0, x=Éb, zeBC, so (3) 
Xx 
dy =— 
Kw) =0, x=0, zeDA, t>0 (14) 
Xx 
E -1 xeCD z=P, t>0 (15) 
Z 


la profundidad del dominio de solución (P) 
debe ser mayor que la máxima posición del 
frente de humedecimiento durante el riego. 

La solución de la ecuación de Richards hace 
indispensable representar las propiedades 
del expresando el 
potencial de presión (y) como una función 
del contenido volumétrico de agua (0) y la 


hidrodinámicas suelo, 


conductividad hidráulica K como una función 


L 


a 


Figura 1. Dominio de solución para la ecuación de Richards bidimensional. 
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de 6. Como es señalado por Fuentes et al. 
(1992), la combinación de las características 
hidrodinámicas de Fujita (1952) y Parlange et 
al. (1982) es conveniente en estudios teóricos, 
como la construcción de soluciones analíticas 
exactas; en estudios experimentales puede 
ser más conveniente la utilización de la 
combinación de la curva de retención propues- 
ta por van Genuchten (1980), considerando 
la restricción de Burdine (1953) con la curva 
de conductividad hidráulica propuesta por 
Brooks y Corey (1964), debido a que satisfacen 
las propiedades integrales de la infiltración 
y a la facilidad para la identificación de sus 
parámetros. 

La curva de retención propuesta por 
van Genuchten (1980) se escribe como: 
[9 (w) - 0/19,- 0] = [1+ (w/w,)"F”, donde 
y, es un valor característico de la presión del 
agua en el suelo, m y n son dos parámetros de 
forma empíricos relacionados por la restricción 
de Burdine (1953): m =1-2/n,con0<m<l y 
n > 2, 0; es el contenido volumétrico de agua a 
saturación efectiva del suelo y 0, es el contenido 
volumétrico de agua residual. La conductividad 
hidráulica propuesta por Brooks y Corey (1964) 
se representa como K(0) = K_[(0 — 0,)/(8, - 0,)], 
donde y es un parámetro de forma empírico y 
positivo. 

La lámina infiltrada se calcula para cada po- 
sición x, donde se resuelven numéricamente 
las ecuaciones de Saint-Venant: 


F 
I(x,,t)=Sl0(x,2,+)-0(x,,2,0)]dz (16) 
0 


donde F es la posición del frente de hume- 
decimiento y 0(x,z,0) es el contenido volumé- 
trico de agua inicial. 


Solución numérica 


El esquema numérico para la fase de avance ha 
sido expuesto por los autores (Saucedo et al., 
2005). La disposición de las celdas de cálculo 
para las fases de avance, almacenamiento, 


consumo y recesión se muestra en la figura 
2. La configuración de las celdas para la fase 
de recesión obedece al hecho de que en lugar 
de utilizar un esquema de paso de tiempo 
constante, se hace uso de un esquema de paso 
de espacio fijo. 

La forma discreta de la ecuación de 
continuidad para las fases de almacenamiento, 
consumo y recesión se escribe como: 


3, +(1- 0)9,J8+- [oq, +(1- (0)q, 18 
7 [9n, +(1- 9), +01, + (1-9), lx, —x1) 
+ [oh,+(1=9)M +01, +10), ](x,—x,)=0 


(17) 


para la ecuación de momentum, la forma 
discreta es: 


l r 
y 
Y 
] m 
: Z r 
ti+1 7 
ti 
] m 
X-1  * Xi¡+1 Xy 


Figura 2. Disposición de las celdas de cálculo para las fases 
de avance, almacenamiento, consumo y recesión. 
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27 H(q,-q/18+ (gh? -72) (Mm, —h ,J8t 
+x,-x,) 09, +(1-0)9,-09, (1 0); 
+25 h*(x,—x,) [0], +(1-0)J, -J,] (18) 
+p3tqh(x,-x)L01,+(1-0)1, 01, -(1-0)1]=0 


donde los coeficientes se calculan a partir de 
los valores pertenecientes al nivel de tiempo 
anterior: 7 = (1 — Q),, + 0q, h = (1 - 0)h,, + 0h, 
(Saucedo et al., 2005). 

Se introducen las variables pequeñas 
(Strelkoff y Katopodes, 1977): 8h,, Sq,, 0h, y 04, 
tales que, para las celdas interiores se tiene: 
h,=h,+06h, h, =h,,+ 0h, q,=q,+ 0q, y q,=4,, + 


9q,. Se sustituyen en las formas discretas de las 


ecuaciones de continuidad y momentum, y se 
obtiene un sistema de ecuaciones algebraicas, 
cuya solución permite avanzar en el tiempo 
los valores de las variables hidráulicas. 


Fase de almacenamiento 


Debido que en la última celda se tiene un 
tirante nulo en la punta de la onda al inicio 
de la fase de almacenamiento, el cual no 
sirve para aproximar el tirante en el tiempo 
siguiente, se emplea solamente en esa etapa, 
denominada de enlace, el tirante que se 
tiene en un punto anterior en el espacio para 
avanzar los perfiles de flujo. La etapa cuasi- 
estacionaria se tiene una vez que el tirante 
en el extremo aguas abajo de la melga es 
diferente de cero, y se emplea para estimar 
el tirante en el tiempo siguiente y avanzar 
los perfiles de flujo. Para la última celda 
en la etapa de enlace se tiene h, = h, + 6h, 
h, = h, + Sh, q, = q, + 9, mientras que para la 
etapa cuasi-estacionaria se escribe h, =h, + 0h, 
h, =h,, + Sh, q,= q, + q; puede apreciarse 
que la forma de las aproximaciones en 
ambas etapas difiere en el cálculo de h,. Las 
expresiones que se obtienen para el sistema 
matricial se muestran en el anexo. 


Fase de consumo o recesión vertical 


La fase de consumo o de recesión vertical 
comienza una vez que se tiene en el dominio 


de solución una cantidad de agua que 
permite disponer de una lámina media igual 
a la lámina de riego que desea aplicarse, y 
termina en el momento en que se consume el 
agua en la cabecera de la melga, dando inicio 
a la fase de recesión. Para modelar la fase de 
consumo se utilizan las aproximaciones de 
la fase de almacenamiento en su etapa cuasi- 
estacionaria, considerando la eliminación 
del gasto en la cabecera de la melga en el 
cálculo de la primera celda. Las entradas del 
sistema matricial para la primera celda se 
muestran en el anexo, mientras que para las 
celdas interiores y la última celda se hace uso 
de las ecuaciones para dichos tipos de celda 
empleadas en la fase de almacenamiento. 


Fase de recesión 


La fase de recesión ha sido modelada 
mediante un esquema con paso de espacio 
constante (0x), que permite evitar los errores 
que se involucran en la reproducción del 
frente de recesión cuando se hace uso de 
un esquema con paso de tiempo constante; 
dichos errores tienen que ver con el hecho de 
que difícilmente se obtiene un tirante nulo en 
una posición de la malla donde se resuelven 
las ecuaciones de Saint-Venant para un 
tiempo que sea múltiplo del paso de tiempo 
constante empleado en la discretización, 
motivando que el frente de recesión se corra 
de manera anticipada al nudo siguiente; no 
obstante que esta aproximación produce 
errores de balance poco significativos si se 
agrega a la lámina infiltrada el tirante que 
permanece sobre el nudo desde el cual se 
efectúa la traslación del frente de recesión, es 
posible mejorar la aproximación utilizando 
un esquema con paso de espacio constante. 
Es necesario realizar una interpolación entre 
las variables tanto de las ecuaciones de Saint- 
Venant como de la ecuación de Richards para 
pasar de la malla de paso de espacio variable 
utilizada en las fases de consumo a una malla 
de paso de espacio constante. 


Saucedo et al., Modelo hidrodinámico completo para riego por melgas 


La utilización de un esquema de paso 
de espacio constante genera un sistema de 
ecuaciones, en el cual el número de incógnitas 
es inferior en dos unidades al número de 
ecuaciones, motivo por el que se acepta no 
resolver la ecuación de conservación de 
momentum en la primera y la última de las 
celdas de cálculo. El paso de tiempo se calcula 
utilizando las variables hidráulicas de la 
primera celda, teniendo en cuenta la forma 
discreta de la ecuación de continuidad: 


St (1 0), = Sxlo1, as ( Y), +1,)-1,) (19) 


es importante señalar que se desconoce a 
priori el valor de 1, motivo por el cual es 
necesario proporcionar un estimador inicial, 
en este caso se tienen en cuenta los valores 
de lámina infiltrada correspondientes al 
nudo de coordenada r en los dos niveles de 
tiempo anteriores, de tal forma que es posible 
aproximar la velocidad de infiltración y calcular 
el primer estimador como 1; =1,, + V,6t, donde 
V, = (15 — 12) /St, siendo k el nivel de tiempo, 
y se resuelve la ecuación (19) para el paso de 
tiempo mediante un procedimiento iterativo. 
Con el valor del paso de tiempo así obtenido 
se avanzan los perfiles de flujo superficial 
y subsuperficial, resolviendo el sistema de 
ecuaciones que se muestra en el anexo. 


Solución numérica de la ecuación de Richards 


La ecuación de Richards bidimensional es 
discretizada en el espacio utilizando ele- 
mentos finitos y en el tiempo mediante un 
esquema implícito en diferencias finitas. El 
procedimiento se encuentra bien documenta- 
do en la literatura y puede consultarse, por 
ejemplo, en Saucedo et al. (2002). La solución 
del sistema de ecuaciones algebraicas que 
resulta de la aplicación del método del 
elemento finito es efectuada usando el método 
del gradiente conjugado precondicionado 
(Noor y Peters, 1987). 


Aplicación 


Con la finalidad de verificar el funcionamien- 
to del modelo desarrollado para las cuatro 
fases del riego por melgas, basado en el 
interno de las 
ecuaciones de Saint-Venant y Richards, se 


acoplamiento numérico 


realiza una aplicación empleado datos de un 
experimento reportado por Fuentes (1992). Se 
dispone de los siguientes datos: gasto unita- 
rio q,= 0.0032 m*/s/m; pendiente topográfica 
J, = 0.002 m/m, longitud de la melga L = 100 
m; parámetros para la ley de resistencia de 
Fuentes et al. (2004) d=1, k=1/54; parámetro 
en la ecuación de cantidad de movimiento 
PB = 2; valor inicial del contenido volumé- 
0.2749 cm*/cm?. Los 
parámetros 0, 0, m y n de las características 


trico de agua 0, = 


hidrodinámicas fueron determinados a partir 
de la curva granulométrica del suelo franco de 
Montecillo. Siguiendo la metodología sugerida 
por Fuentes (1992) resulta 0, = 0.4865 cm? /cm, 
9, = 0.0 cm*/cm?, m = 0.1258, n = 2.2878 y 
n = 11.0. Los valores de los parámetros 
y, =-32.75 cm y K_=1.84 cm/h se obtuvieron 
mediante la aplicación de un procedimiento 
inverso para describir datos de una prueba 
de 
mencionado. 


riego efectuada en el experimento 

Se ha utilizado un paso de tiempo 
constante df = 0.5 s tanto para las ecuaciones 
de Saint-Venant como para la ecuación de 
Richards. Para la ecuación de Richards se 
utilizaron pasos de espacio mínimo y máximo: 
Az, =0.02 cm y Az, . =1.0cm, Ax, . =0.02 cm 
y Ax 
guarda semejanza con las reportadas en la 
literatura, e.g., Strelkoff y Katopodes (1977): 
St. = 5 Ss, Akanbi y Katopodes (1988): 


mín 


= 1.0 cm. La discretización utilizada 


máx 


db =15, para las cuales se reportan resultados 
consistentes en las variables hidráulicas del 
riego. La discretización para la ecuación de 
Richards es más fina que la reportada por 
Saucedo et al. (2002), quienes comparan perfiles 
de humedad proporcionados por una solución 
analítica exacta con los obtenidos mediante 
una solución numérica para la ecuación de 
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Richards; los valores del error relativo máximo 
para el caso de dos suelos contrastantes: una 
arcilla clara de Yolo, Estados Unidos, y una 
arena del río Isére, Francia, presentan una cota 
superior de 0.25% cuando se utiliza un paso 
de tiempo máximo St. = 30 s y pasos de 
espacio mínimo y máximo Az, = 0.04 cm y 
Az... = 10.0 cm. El suelo franco de Montecillo 
tiene un comportamiento hidrodinámico in- 
termedio entre los que presentan los suelos 
tratados por estos autores; en consecuencia, 
con la discretización utilizada, el error relativo 
máximo debe tener como una cota superior 
0.25%. 

En la figura 3 se muestran los perfiles de 
flujo obtenidos durante la fase de avance y 
en la figura 4, las curvas de avance y recesión 
para la melga cerrada. En la figura 5 se tiene 
un detalle del comportamiento de los tirantes 
de agua en la parte final de la melga durante 
la fase de almacenamiento. La figura 6 
muestra las cuatro fases del riego por melgas. 
En la figura 7 se representa la evolución y 
recesión vertical del tirante en la cabecera de 
la melga vista a través del comportamiento 
del potencial de presión. En la figura 8 se 
tiene la evolución del potencial de presión 
proporcionado por la ecuación de Richards 
bidimensional a lo largo de la melga durante 
la fase de recesión. 

En la figura 9 se muestra el buen ajuste 
que presenta la modelo 
hidrodinámico a datos de avance registrados 


respuesta del 


en el experimento reportado por Fuentes 
(1992). 


Conclusiones 


Se ha desarrollado un modelo numérico para 
riego por melgas basado en el acoplamiento 
interno de las ecuaciones de Saint-Venant y 
Richards, el cual permite describir las fases 
de avance, almacenamiento, consumo y 
recesión del riego por melgas. La solución de 
las ecuaciones de Saint-Venant se aproxima 
utilizando un esquema Lagrangiano en 


diferencias finitas, mientras que para la 


ecuación de Richards se utilizan elementos 
finitos para la integración en el espacio 
finitas 


y diferencias implícitas 


integración en el tiempo. Los resultados 


para la 


obtenidos permiten concluir que el modelo 
puede ser utilizado con fines de diseño de 
riego por melgas. 


Anexo 
Fase de almacenamiento 


Las entradas del sistema matricial relacionadas 
con la solución de la forma discreta de la 
ecuación de continuidad para el caso de la 
primera celda son: 


lo(x,-x)151,+[008:]59,+[1-0)x,-x,)]5h, 
(A.1) 
1)+41-9)14, 1) (x,-x,) 


=5Hg; q) +90, 


para las celdas interiores se tiene: 


Eostlag,+ lo(x,- x,)] 3h, + [08:]57, 


+[41-oJx,—x]5%,=8+(q¡-q,)  (A2,) 


+[9(,-1,) +10, 1, x,-x,) 


y para la última celda, en la etapa de enlace: 


[vst]5q,+ Eo(x, — 1181, +1 41-0)x,-x,)]8h, 


= —5t9,-L[9(1,-1,)+(1-0)01,+1,)](x,-x,) (A.3) 
y en la etapa cuasi-estacionaria: 


[ost ]3g,+Ep(x, -1]00,+-41=01x,-x,)]5, 


= —5tq, 190, -1,)+(1-90, +1)1(,-x) (4.4) 
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Tirante /lámina infiltrada (cm) 


10 (00 


Distancia (cm) 


Perfiles de flujo obtenidos durante la fase de avance, correspondientes a tiempos de 1, 5, 10 y 20 minutos. 
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Figura 4. Curvas de avance y recesión. 
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Distancia (cm) 


Figura 6. Aspectos de las cuatro fases del riego por melgas. 
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Potencial de presión (cm) 
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Figura 7. Evolución del potencial de presión en la cabecera de la melga. 
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Figura 8. Evolución del potencial de presión a lo largo de la melga (fase de recesión). 
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Figura 9. Ajuste a datos de avance de una prueba de riego. 


Las entradas del sistema matricial relacio- 
nadas con la solución de la forma discreta de 
la ecuación de momentum para el caso de la 
primera celda son: 


coso 7?) - 3088 (x,- pon 


] 


OS A 


mM, 


+ |w0t (gh q>)-3(1- d)e Sth * (x,— ste sh, 


m 


=2qh8t (q, -q,,) +54 gh? 7), hi) 
got (x,-x,) 1-0), +0) —J,] 


-Bqh(x,-x,)19,-1)+01-910,-1,)] (45) 


para las celdas interiores se tiene: 


- 27 host+4h (x,-—x,)+ 7 Mx, -x,) ] a, 
q; 


.- ost gh*- 77) -3035H x, 00h + 
J 


gio (xx) + (1 re, an, 


ES [og 7) 3(1-0)2 5H (x,-x,) Lo, 
= 2qh0t(q, = Am) +5igh E 7) (h; E h,,) 
- gm (x, —x 1003), +0 J,=J,] 


-pgix,- 1) [9(,-1)+(0-0)4, —1,)) (A.6) 
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Para la última celda, en la etapa de enlace se 
escribe: 


- 27 ho8t+01 Ax, ro (x,- 0 a 


+ as(sr-5)- 303 5th"(x, 90d ln 
11 (A7) 


+ lw0sH( gh? 72)]6n, =27H (1-0)8tg; 
+0t(1— gh? - q),- hA1-0NMx, — xq; 
25H, 100; -J,) BGM, -x)1d 0)1,+00,-1)] 


y en la etapa cuasi-estacionaria: 


ños t+ ph, xp) + AS xy) La 
q 


+ -oi(gi-47) -3 q Sth Ax, =p lr, 
h; 


+[wst(gn?-32]5n, =-27h(1- 0)8tg,, 


+8t(1- og 7)M,-h,) 35 (x,- xp) 1), -J,) 
- pq h(x,- xl -0)1,+0(1,-1,) (A.8) 


Fase de consumo 


Las entradas del sistema matricial relacionadas 
con la solución de la forma discreta de la 
ecuación de continuidad para la primera celda 
son: 


lb(x, —x,) ]5n, + lost l5g,+[(1-91x,-x,) 15h, 


=— 8tq y +L0(, -1,)+ (1-0) 1,96%, — x,) (A.9) 


Las entradas del sistema matricial relacio- 
nadas con la solución de la forma discreta de la 
ecuación de momentum para la primera celda 
son: 


E (en? - q) 3035th*(x, — sol 
j 


.[a0 5070, 100701, 


+ ds 3-72) -3(1-0)25tHh (x,—x,) La 
=-23h8tq , + Skgh? 3?) (1h) (A.10) 
- goth *(x,- xa > IA +0 J; Ji) 

-Bah(x, —x)Í9(, 1) + (1-9), —1,)] 

Fase de recesión 

Las entradas del sistema matricial relacionadas 
con la solución de la forma discreta de la 


ecuación de continuidad para la segunda celda 
son: 


[051] 57,+ [(1- 0)5x ]5,=-5tg,, 


(A.11) 
+[9(,-1,)+(1-0)0,,-1,)]6x 


Para las celdas interiores se tiene: 


E ost] 59,+[v5x ]51,+ [ost]5g, + [(1-0)5:]5n, 


=54(q,—qn)+L0(,-1,)+(1= 910, -1,)]óx (A.12) 
Para la última celda se escribe: 


[os+J5g,+[Egár 151, + [-(1-0)5x Jn, =-8tq, 


[94,-1,)+1-0)4,-1)]65x A.13) 


Las entradas del sistema matricial relacio- 
nadas con la solución de la forma discreta de 
la ecuación de momentum para la segunda 
celda: 
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Ey 5t+(1-0)/178x +9 3 E Si 


m 


+ Loss 3-72)-3(1-9)g8t8xT? la] Sh, 


m 


= 27 h6t (q; Sl 4m) + sH(gh- 77) (h; 5 h) (4.14) 


—gotóxli? (1-0) J,, + 03, -J,) 
= Bqhóx [o(1, —1,) H= 9, Lg )] 


Para las celdas interiores se tiene: 


-2q ho8t+ 01 *8x +4 a l 59, 
4 


+|- 08H gh?- 77) 3090 Sh, 
j 


+ agro (1-0)? 8x + a e 57, 


m 


+ alg - q) 3(1-0)25t8xh* le sh, 


m 


e 27 h8t(q, — Gm ) + st(gn*-7)m, A) 
got? [1—0)] + 0J,-J, ) 


-parixlod, 1) +19), —1,,)] (A.15) 


Recibido: 17/02/09 
Aprobado: 16/08/10 
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Abstract 


SAUCEDO, H., ZAVALA, M. és FUENTES, C. Complete hydrodynamic model for border 
irrigation. Water Technology and Sciences, formerly Hydraulic engineering in Mexico 
(in Spanish). Vol. II, No. 2, April-June, 2011, pp. 23-38. 


A numerical model for border irrigation is presented —based on the internal coupling of Saint- 
Venant and Richards equations —that makes it possible to describe the advance, wetting, 
depletion, and recession phases. The solution of Saint-Venant equations is approached by a 
finite-difference Lagrangian scheme, while for Richards equation finite elements for the space 
discretization and implicit finite-differences for the time discretization are used. The results 
obtained with the model application lead to conclude that the model can be used for border 
irrigation design. 


Keywords: advance phase, wetting phase, depletion phase, recession phase, internal numerical 
coupling. 
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